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Abstract. We investigate the effect of a rotation of the In- 
terplanetary Magnetic Field (IMF) on the transport of mag- 
netospherie ion populations at Mercury. We focus on ions 
of planetary origin and investigate their large-scale circula- 
tion using three-dimensional single-particle simulations. We 
show that a nonzero Bx component of the IMF leads to a 
pronounced asymmetry in the overall circulation pattern. In 
particular, we demonstrate that the centrifugal acceleration 
due to curvature of the E x B drift paths is more pronounced 
in one hemisphere than the other, leading to filling of the 
magnetospheric lobes and plasma sheet with more or less 
energetic material depending upon the hemisphere of origin. 
Using a time-varying electric and magnetic field model, we 
investigate the response of ions to rapid (a few tens of sec- 
onds) re -oriental ion of the IMF* We show that, for ions with 
gyroperiods comparable to the field variation time scale, the 
inductive electric field should lead to significant nonadiabatic 
energization, up to several hundreds of eVs or a few keVs. It 
thus appears that IMF turning at Mercury should lead to lo- 
calized loading of the magnetosphere with energetic material 
of planetary origin (e.g„ Na + )« 

Keywords, Magnetospheric physics (Planetary magneto- 
spheres) 


1 Introduction 

With an orbit comprised between M).3 and ^0.47 astronomi- 
cal units, Mercury is located in a region of the Parker’s spiral 
where the radial component of the Interplanetary Magnetic 
Field (IMF) is dominant. Because of its appreciable though 
weak (about 1000 times smaller than that of Earth) intrin- 
sic magnetic field, a miniature magnetosphere forms around 
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the planet (e.g., Ness, 1979; Russell and Walker, 1985; An- 
derson et al., 2008) that is expected to be very dynamical 
and subjected to prominent reconfigurations due to variations 
of the IMF or variations of the solar wind dynamical pres- 
sure (e.g., Slavin et al., 2010), In a recent study (Delcourt 
et al., 2010), we examined the response of magnetospheric 
ions to compression of the magnetosphere due to abrupt so- 
lar wind pressure increase. In that study, we focused on 
ions originating from the planetary exosphere that are sus- 
pected to significantly contribute to the inner hermean mag- 
netotail as evidenced from MESSENGER data (Zurbuchen 
et aL, 2008). We showed that, during such events, a mass 
selective nonadiabatic heating is at work due to the inductive 
electric field. In the present study, we investigate the effect of 
IMF re-orientation on the dynamics of magnetospheric ions 
and show that a similar heating should operate. In Sect, 2, 
we describe the procedure adopted to simulate the ion trans- 
port during turning of the IMF. In Sect* 3, we discuss some 
results of trajectory computations in steady state, while those 
obtained during dynamical reconfigurations of the magneto- 
sphere are presented in Sect. 4. 


2 Modeling of magnetosphere reconfiguration 

To examine the response of magnetospheric ions to IMF 
turning at Mercury, we performed single-particle simula- 
tions in three-dimensional time- varying electric and mag- 
netic fields that provide a first-order description of the lier- 
mean environment. As for the magnetic field, it was ac- 
counted for using a modified version of the model of Luh- 
mann and Friesen (1979) (hereinafter referred to as LF-79)* 
This simple model that reproduces basic features of the mag- 
netosphere, consists of a planetary dipole field and an in- 
finite Harris-type current sheet (see Delcourt et al., 2010). 
In the present study, we take for the planetary dipole mo- 
ment M = 250 nT /?^ (e.g., Anderson et al., 2008). We also 
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Fig. 1* Magnetic field lines in the noon-midnight meridian plane 
for IMF Bx = 0 and IMF Bx = 20 nT (color-coded in blue and red, 
respectively). The latitude of the magnetic field line foot on the 
planet varies by step of 5°. The X-axis points toward the tail, the 
Y-axis, from dusk to dawn. 

assume a current sheet field, Bf* and half- thickness, L y of 
50 nT and 0.5 Rm, respectively, which leads to a subsolar 
magnetopause near 1.5 Rm. On the other hand, we simulate 
the effect of the IMF by superposing a Bx component to the 
above modified LF-79. Specifically, we consider a situation 
where the IMF Bx varies from 0 up to 20 nT, The magneto- 
spheric configurations obtained in the initial and final states 
are shown in Fig, i that presents magnetic field lines in the 
noon-midnight meridian plane (in this figure as well as in the 
following, note that the X-axis points in the aniisunward di- 
rection, and the Y-axis, from dusk to dawn). The warping of 
the magnetosphere that results from a non- zero IMF Bx is 
clearly apparent from Fig. 1 , 

Like in previous studies (e.g., Del court et ah, 2002, 2010), 
the large-scale magneto spheric convection at Mercury' is 
modeled by assuming a situation similar to that at Earth, with 
a two-cell pattern of potential distribution in the immediate 
vicinity of the planet as featured in the Volland-Stem descrip- 
tion (see, e.g,, Volland, 1978). The inclusion of a nonzero 
IMF Bx however leads us to implement some modification. 
Indeed, in the initial state with IMF Bx = 0 (blue lines in 
Fig. 1), the symmetry between Northern and Southern Hemi- 
spheres allows us to use the same potential distribution in 
both hemispheres. This is not the case in the final state with 
IMF fix“2GnT (red lines in Fig. 1), In this latter case, dis- 
tinct potential distributions must be considered which would 
otherwise lead to unrealistically abrupt variations of the elec- 
tric field (hence, unrealistic variations of the particle mag- 
netic moment) at the transition between one hemisphere and 


NORTH - IMF S, = 0 NORTH - IMF B, = 20 nT 




SOUTH - IMF B* = 0 SOUTH - IMF B, = 20 nT 




Fig, 2. Potential distribution in (top) Northern and (bottom) South- 
ern Hemispheres for (left) IMF B\ — 0 and (right) IMF Bx — 20 nT. 
Circles indicate latitudes by steps of 10°. Equi potentials are shown 
by steps of 2kV> 

the other. The potential distribution in the Volland (1978) 
model is characterized by two parameters, viz., (i) the co- 
latitude of the Convection Reversal Boundary (CRB) above 
and below which plasma is flowing in antisunward and sun- 
ward directions, respectively, (ii) the cross-polar cap poten- 
tial drop. As for the CRB colatitude, it is evaluated at each 
step of the simulation by calculating in each hemisphere the 
latitude of the last closed field line in the noon meridian 
plane. For the case shown in Fig. I (red lines), this colati- 
tude is of the order of 58° and 43° in Northern and Southern 
Hemispheres, respectively. As for the cross-polar cap po- 
tential drop, for a given value in the Northern Hemisphere, 
the value in the opposite hemisphere is calculated in such a 
manner that the E x B drift speeds obtained in the distant 
magnetotail from one and the other potential distributions 
are identical. The resulting distributions can be appreciated 
in Fig. 2 that shows equipotenlials in both hemispheres (top 
and bottom rows) in the initial and final states (left and right 
panels). While the left panels (IMF Bx = 0) of Fig, 2 display 
similar patterns, notable differences can be seen in the right 
panels for IMF Bx — 20 nT, with tailward flow over a wider 
latitudinal range in the Northern Hemisphere (top) than in the 
southern one (bottom). 

Figure 2 relates to the convection electric field that results 
from coupling with the solar wind. On the other hand, the 
transition from IMF Bx = 0 to IMF Bx = 20 nT gives rise 
to an inductive electric field that must also be accounted for 
This very' transition was modeled in a like manner to Dd- 
court et ah (2010). At a given position r and a given time i 
of the transition, the magnetic field was obtained as: 
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Fig. 3. (Left) Evolution of a given magnetic field line from IMF 
£?X — 0 (blue color code) to IMF fix = 20 nT (red color code). In 
both hemispheres, the black line shows the E x B drift path of a 
zero energy particle during the IMF turning. (Right top) Electric 
field and (right bottom) E x B drift speed versus time. The time 
scale of the IMF turning is set to 20 s. 

B(r,t) = Bo(rnmiBf(r)-BUr)] ( 1 ) 

where B 0 and B / correspond to initial and final configu- 
rations, respectively. In Eq. (!), /(f) describes the rate of 
change of the magnetic field. As described in Delcourt et 
ah (2010), imposing a zero induced electric field as well as 
its zero slope at the onset and at the end of the magnetic tran- 
sition leads to a polynomial of degree 5 for /(f). Denoting 
by tb the time scale of the magnetic transition, this polyno- 
mial smoothly varies between 0 at t = 0 and 1 at t ^ Tg . 

As for the time-dependent electric field, it w r as consid- 
ered to be the sum of two contributions. The first contribu- 
tion is that induced by the time varying magnetic field, vi 2 ., 
E = —BA/Bt where the vector potential A is such that curl 
( A)= B as given by Eq. (1). Since the vector potential is un- 
determined to the gradient of a scalar field, a second contri- 
bution is also considered that cancels the parallel compo- 
nent of the induced electric field, that is, (0A/9r) B = 0. The 
calculation of this second contribution is described in Ap- 
pendix B of Delcourt et at, (1990). As discussed by Heikkila 
and Pellinen (1977), this latter contribution actually redis- 
tributes the electric field in the perpendicular direction and 
accounts for plasma polarization. At a given position r and 
time t , the electric field due to IMF turning is thus obtained 
as: 

( 2 ) 

where V$poi is the electric potential due to the free charges, 
and A is the instantaneous vector potential such that: 

A(rj) = A,,(r) + /(0[A/(r) - A*(r)] (3) 

where A 0 and A / relate to initial and final configurations, 
respectively. 

The electric and magnetic field variations obtained from 
Eqs. (1-3) can be appreciated in Fig. 3 that shows the evolu- 
tion of two field lines anchored at the same latitude (75°) in 


the Northern and Southern Hemispheres, In the left panel of 
Fig. 3, it can be seen that, as a result of the IMF fix growth, 
the field line in the Northern Hemisphere is rapidly brought 
sunward while that in the Southern Hemisphere is pulled tail- 
ward. Note that, in this figure as well as in the following sim- 
ulations, the time scale tb of the magnetic transition was set 
to 20 s. Because propagation effects are not accounted for in 
the present model, this 20 s time scale describes the duration 
of the local magnetic field reconfiguration, the overall recon- 
figuration of the magnetosphere itself occurring on a larger 
time scale (a few minutes). To further illustrate the electric 
field variation during this magnetic transition, the black lines 
in the left panel of Fig. 3 show the trajectories of zero en- 
ergy particles that experience only E x B drift from initial 
to final configurations. The right panels of Fig. 3 show the 
electric field and E x B drift speed along these paths. Con- 
sistently with the field line evolution displayed in the left 
panel, a dawnward induced electric field is obtained in the 
Northern Hemisphere (solid line) with a peak magnitude of 
—3 mV m” 1 ; hence, a peak E x B drift speed of —50 km s“ 1 
(average value of the order of 30kms _1 ). In the Southern 
Hemisphere, a duskward induced electric field is obtained 
that leads to a peak E x B drift speed somewhat smaller than 
in the Northern Hemisphere, of the order of 35 km s” 1 . 

Once electric and magnetic fields were defined, test ions 
were traced throughout Mercury's magnetosphere using the 
full equation of motion (integrated by means of a fourth- 
order Runge-Kutta). Such a tracing technique is fairly stan- 
dard and has been used for instance in a variety of studies of 
plasma transport in the terrestrial magnetosphere. In order 
to examine the impact of changes in the solar wind on the 
large-scale magnetospheric circulation, some of these stud- 
ies do not use analytical or semi-empirical fields but results 
of MHD simulations (e.g., Moore et ak, 2005; Peroomian 
and El-Alaoui, 2008), Such an approach is also used for spe- 
cific case studies (e,g, ? Ashour-Abdalla et aL, 2010). In this 
latter approach, the electric and magnetic fields at the test 
particle position are interpolated from the three-dimensional 
grid of the MHD simulations. The aim of the present study 
is not to examine specific solar wind cases but to investigate 
from a general perspective some transport features associ- 
ated with IMF rotation; hence, the use of a simple analytical 
description for both electric and magnetic fields. It should 
he noted that, be they analytical or MHD-based, both ap- 
proaches are not self-consistent in the sense that there is no 
iteration and feedback between the calculated particle mo- 
tion and the field. Such approaches are nonetheless use- 
ful to understand the complexity of magnetospheric parti- 
cle dynamics. In this regard, the three-dimensional electric 
and magnetic fields used In the present study reproduce ba- 
sic magnetospheric features expected at Mercury, with anti- 
sunward transport from the cusp region over the polar cap 
and return flow toward the sun at low latitudes. Note also 
that the present trajectory computations do not include small- 
scale acceleration processes. The energization affecting ions 
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Fig. 4. Model Na + trajectories for (top) IMF fix =0 and (bottom) 
IMF fix = 20 nT. (from left to right) Trajectory projections in the 
X-Z planes, energy, and magnetic moment (normalized to the ini- 
tial value) as a function of time. The ions are launched with lOcV 
energy and 30° pitch angle on the noon meridian and from different 
latitudes (from 60° to 75° by steps of 5°) in the Northern Hemi- 
sphere (blue color code) and in the southern one (ydtow-red color 
code), 

during their three-dimensional transport is solely due to the 
large-scale convection electric field as well as to the surging 
electric field associated with IMF rotation (the role of grav- 
ity being negligible at Mercury). The modeling technique 
adopted here (Eqs. 1-3) ensures that the parallel electric field 
vanishes regardless of the particle position in the magneto- 
sphere, This parallel electric field is not arbitrarily set to 
zero. As described above, the condition E B ~ 0 is achieved 
by assuming that, in first approximation, magnetic field lines 
are equipotential and by making use of a polarization elec- 
tric field, VOpoj, that redistributes the field in the perpendic- 
ular direction {see Appendix B of Delcourt et al. (1 990) for 
a detailed description of the modeling technique). As will be 
shown in Sect. 4 hereinafter, explicit calculation of the paral- 
lel component of the three-dimensional electric field demon- 
strates that it is indeed negligible (less than 10 -4 mV m -1 
along the computed ion trajectories). 

3 Centrifugal acceleration in magnetospheric lobes 

In the adiabatic (guiding center) approximation, the parallel 
equation of motion in the absence of parallel electric field 
and neglecting the gravitational acceleration writes (see, e.g., 
equation 1.20 of Northrop, 1963): 

Bb 0b . a 0B 

^ = ^. [ - + V // - + VVV* ] -£- (4) 

Here, V g is the E x B drift velocity, V//* the parallel speed 
of the guiding center, b , a unit vector in the B direction, s , 
the curvilinear abscissa along the field line, and m and ju, 


Fig. 5, Identical to Fig, 4 but for protons. 


the particle mass and magnetic moment, respectively. Equa- 
tion (4) displays two distinct curvature related acceleration 
terms, viz,, the second term on the right hand side of Eq. (4) 
is due to curvature of the magnetic field lines which can lead 
to prominent acceleration upon crossing of a field reversal as 
is the case in the magnetotaii current sheet. The third term 
on the right hand side of Eq. (4) is a centrifugal acceleration 
that depends upon curvature of the E x B drift paths. In the 
Earth’s magnetosphere, this centrifugal effect can lead to up- 
ward ton energization of a few eV$ or a few tens of eVs at 
low altitudes (e.g., Cladis, 1986) or to focusing (referred to 
as ''centrifugal trapping") of slow plasma sheet particles into 
the equatorial region (e.g,, Delcourt ef a)., 1993). 

It was shown in Delcourt et al. (2002) that, because curva- 
ture radii of the E x B drift paths are much smaller at Mer- 
cury than at Earth, the above centrifugal effect (third term in 
Eq. 4) is enhanced and may lead to quite substantial paral- 
lel energization during transport from high to low latitudes. 
As a matter of fact, it was demonstrated that, as a result of 
this effect, heavy ions originating from Mercury's exosphere 
(e.g., Na + ) may experience energization of several tens or 
a few hundreds of eV within minutes, suggesting a loading 
of the magnetospheric lobes more rapid and more energetic 
than at Earth. This efficient centrifugal energization of low- 
energy planetary material can be appreciated in Figs. 4 and 5 
that present results of trajectory computations for both Na + 
and IP, respectively. In these figures, test ions are launched 
from different latitudes (from 60° to 75° by steps of 5°) in 
both hemispheres in the dayside sector. 

The top rows of Figs. 4 and 5 show the results obtained 
for IMF fix = 0. In this case, symmetrical patterns are ob- 
tained regardless of the hemisphere of origin. Upon transport 
into the lobes, Na + ions experience a net energy gain of sev- 
eral hundreds of e V (top center panel of Fig. 4) while that of 
protons is of the order of 60 eV (top center panel of Fig. 5). 
Subsequently, the ions intercept the equatorial neutral sheet 
where their motion turns non adiabatic as evidenced from the 
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Fig* 6. (top) Model normalized density and (bottom) average energy 
of Na + ions in the noon-midnight meridian plane for (left) IMF 
B x = 0 and (right) IMF B x = 20 nT. Density and energy are coded 
according lo the color scales at right. 


targe pi changes in the rightmost panels* During this nona- 
diabatic sequence that is due to the small spatial scale of the 
magnetic field as compared to the ion Larmor radius (e*g„ 
Spetser* 1965; Chen and Palmadesso, 1986: Biichner and Ze- 
lenyi, 1989), Na* rapidly travel westward and are eventually 
lost into the dusk magnetopause. As for protons, they remain 
temporarily trapped at low latitudes before intercepting the 
dusk magnetopause. For both Na^ and protons, it can be 
seen that this rapid westward drift toward the dusk magne- 
topause goes together w ith prominent energization, as appar- 
ent from the abrupt energy increase near the end of the trajec- 
tories in the center panels of Figs, 4 and 5, During oscillation 
about the equator and simultaneous dusk ward transport, the 
net energy gain realized via crossing of the convection elec- 
tric field equipotentials is of the order of 10 kV, i*e*, half of 
the cross- tail potential drop. 

Because it introduces a hemispherical asymmetry, we ex- 
pect significant deviations from these results in the case of 
a nonzero IMF B X . This is illustrated in the bottom rows of 
Figs. 4 and 5 for Na + and H + , respectively* In this latter case, 
it can be seen that, regardless of mass, the test ions originat- 
ing from the southern high- latitude day side sector (trajecto- 
ries coded in yellow-red) travel somewhat further out into 
the magnetotail. Most notably, the E x B related centrifugal 
energization realized along these latter paths is significantly 
lesser than that of ions originating from the Northern Hemi- 
sphere (trajectories coded in blue), of the order of a few hun- 
dreds of eV instead of M keV for Na + (Fig. 4), and a few 
tens of eVs instead of — lOOeV for protons (Fig. 5), The 
net time of flight from the planet surface until the nightside 


Fig, 7. Identical to Fig. 6 but for protons. 

neutral sheet is accordingly larger for Southern Hemisphere 
originating ions. 

A more general view of these different centrifugal effects 
can be obtained from Figs* 6 and 7 that show the results of 
systematic trajectory computations of Na + and protons orig- 
inating from the high-latitude dayside exosphere. The top 
and bottom panels in these figures show the normalized den- 
sity and average energy, respectively, for (left) IMF Bx = 0 
and (right) IMF Bx = 20 nT These figures clearly illustrate 
the E x B ion transport all the way from the dayside sector 
into the equatorial plasma sheet in the nightside sector over 
the polar cap* Such a large-scale transport leads to a sub- 
stantial filling of the magnetospheric lobes prior to panicle 
interaction with the tail current sheet. Note the gradual de- 
crease of the lobe density along these paths, by about 2 orders 
of magnitude as compared to that at low altitudes (top pan- 
els). Not surprisingly, symmetrical patterns are obtained in 
the case of IMF Bx = 0, with lobe energization up to M keV 
for Na^ and MQOeV for protons. As expected from Figs. 4 
and 5, a nonzero IMF Bx leads to a quite distinct pattern with 
prominent differences between northern and southern lobes. 
In this latter case, regardless of ion mass (Figs. 6 and 7), the 
northern lobe exhibits significantly lesser density and larger 
average energy (equivalently, smaller loading time) than the 
southern one, by about one order of magnitude. Also, be- 
cause protons originating from the Southern Hemisphere in- 
tercept the equator further out in the magnetotail, a lesser 
filling of the innermost magnetosphere is obtained. Fi- 
nally, regardless of the IMF orientation, note the prominent 
energization (up to several keVs) realized upon nonadiabatic 
interaction of the ions with the tail neutral sheet ; hence, the 
build-up of a layer of energetic material. 
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Fig* 8* Model Na + trajectories during a 20-s turning of the IMF from Bx = 0 ro fix — 20 nT The test ions are launched from the same 
location (60° latitude on the noon meridian) with 10 eV energy and 30* pitch angle. Top and bottom panels relate to ions originating from 
Northern and Southern Hemispheres, respectively. (From left to right) Trajectory' projections in the X-Z plane, energy, and magnetic moment 
(normalized to the initial value) as a function of time {measured from the onset of the IMF turning). The different colors correspond to 
distinct times of ejection (by steps of 1 min). Color-coded closed circles in the left panels show the position of the particles at the onset of 
the IMF turning* Vertical dotted lines in the right panels delimit the 20-s magnetic transition. 


4 Ion dynamics during IMF turning 

The results presented in Sect. 3 relate to steady IMF orien- 
tations* To illustrate the behavior of ions during turning of 
the IMF, Figs. S and 10 show test ton trajectories using the 
modeling technique described in Sect. 2, with a gradual tran- 
sition (over a time scale of 20 s) from IMF Bx = 0 to IMF 
= 20 nT* In these figures, Na + and protons were launched 
from the same location in the high-latitude dayside sector at 
distinct (col or- coded) times (by steps of 1 min) before the 
onset of the magnetic transition, in the left panel of Figs. 8 
and 10, this onset is indicated by (color-coded) closed cir- 
cles. The top row in these figures show ions launched from 
the Northern Hemisphere, while the bottom row relates to the 
southern one. 

Looking first at Na* in Fig. 8, it can be seen in the left 
panels that, after convection over the polar cap into the night- 
side sector, the test ions reach the equatorial magnetotail 
where they experience non adiabatic interaction with the neu- 
tral sheet* In this figure, the trajectories coded in blue cor- 
respond to ions reaching the equator before the onset of the 
IMF turning. In this case, similar results are obtained regard- 
less of the hemisphere of origin, be it for the trajectory apex 
or for the net fi change. For ions leaving the exosphere at 
somewhat later times, substantial differences are noticeable. 
As discussed in Figs* 4 and 5, it can be seen that as a result of 
the growing IMF Bx, the trajectory' apex of Southern Hemi- 
sphere originating ions gradually moves tailward. For these 
ions, lesser centrifugal accelerations also are achieved in the 
magnetospheric lobes. Most notably, in the right panels of 
Fig* 8, an increase of magnetic moment and concurrent in- 
crease of energy can be seen during the magnetic transition. 



TIME (min) 

Fig* 9* Electric field variations along the trajectories of Na + 
launched 120 s before The onset of IMF turning (trajectories coded 
in yellow in Fig, 8): (top) electric field component in the Y- 
direction, (bottom) parallel electric field (in absolute value) versus 
time. Blue and red colors correspond to Na" launched form North- 
ern and Southern Hemispheres, respectively. Vertical dotted lines 
delimit the 20-s magnetic transition* 

The fj, increase is most pronounced (by two orders of magni- 
tude) for the Na + located at mid-altitudes in the northern lobe 
at the onset of the magnetic transition (trajectories coded in 
green and yellow). This behavior contrasts with that obtained 
for ions located at low altitude (trajectory coded in red) or in 
the equatorial region. 
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Fig. 10, Identical to Fig. 8 but for proions. 


The electric field variations experienced by the ions can 
be appreciated in Fig. 9 that displays the component in the 
dawn- dusk direction throughout transport (top panel). The 
bottom panel of Fig. 9 also shows the magnitude of the par- 
allel component obtained in these computations. The elec- 
tric field profiles presented here are those obtained for the 
test Na + launched two minutes before onset of the IMF turn- 
ing (i.e., trajectories coded in yellow in Fig. 8). Looking at 
the top panel of Fig. 9, it can be seen that, since one initially 
has IMF Z?x = 0* identical Ey variations are first obtained for 
Northern (coded in blue) and Southern Hemisphere (coded in 
red) originating ions, with a magnitude up to at 

low altitudes. From / = 0 until t = 20 s T the IMF Z?x increases 
up to 20 nT and rapid Ey variations of a few mV m _ 1 can 
be seen. These variations lead to a net decrease of the con- 
vection electric field in the Northern Hemisphere and a net 
increase in the southern one, consistently with the magnetic 
field line motion portrayed in Figs. 1 and 3. Subsequently, 
the Ey profiles obtained for the two test Na 4 ‘ are different 
due to their distinct trajectories (see Fig. 8). On the other 
hand, the bottom panel of Fig. 9 confirms that the modeling 
approach adopted in the present study leads to a negligible 
parallel electric field component throughout transport, with 
a magnitude in the 10^ 5 -10 -4 mVm" 1 range (equivalently, 
leading to a net change of parallel energy of at most a few 
Volts during transport). 

A result similar to that displayed in Fig. 8 was reported in 
the study of Delcourt et al. (2010) that was dedicated to the 
dynamics of tons during compression of Mercury’s magne- 
tosphere (see, e.g., Fig. 3 of that study). It was shown in this 
latter study that the above nonadiabatic behavior results from 
resonance between the ion gyro motion and the time- varying 
magnetic field. That is, for ions with zc <$£ z$ (tc being the 
ion cyclotron period), the magnetic field does not change sig- 
nificantly during a cyclotron turn and the motion is adiabatic 
(/t conserving) (e.g. ? Northrop, 1963), In Fig, 8, this is the 


case at low altitudes where the magnetic field is large. On 
the other hand, for ions with r c ^ the work done by the 
induced electric field during a fraction of cyclotron turn is 
not significant enough to yield substantial perpendicular en- 
ergization, a situation that is achieved for ions located in the 
equatorial region at onset. In between (i.e., at mid- latitudes), 
for ions with tc ^ vr , nonadiabatic behavior with possibly 
large fi enhancement is obtained due to phasing between the 
gyromotion and the surging electric field. This nonadiabatic 
energization is controlled by the magnitude of the peak in- 
duced electric field. It thus depends upon r b as well as A A 
(Eq. 3), a shorter time scale or a wider magnetic transition 
resulting into a larger energy gain. 

While the study of Delcourt et al. (2010) investigated the 
effect of a sudden compression of the magnetosphere (mod- 
eled via an abrupt increase of the magnetospheric lobe field, 
the present study is focusing on the impact of IMF turn- 
ing and consequent warping of the magnetosphere. It was 
shown in the above Figs. 4 and 5 that this latter process leads 
to hemispherical asymmetry. In this regard, it can be seen by 
comparison of top and bottom panels in Fig, 8 that the ion 
behaviors do significantly differ depending upon whether the 
particles originate from Northern or Southern Hemispheres. 
During the IMF change considered here (from Bx — 0 to 
Bx — 20 nT), ions exhibit weaker nonadiabatic features in the 
Southern Hemisphere as compared to the northern one. Such 
an impulsive nonadiabatic transport at high latitudes likely 
Forms a specific feature of the hermean environment because 
of both small field variation time scales (small r^) and weak 
planetary magnetic field (large re). 

Since the nonadiabatic response of the panicles depends 
upon gyration period-to-field variation time scale ratio, one 
expects the dynamics of protons to significantly differ from 
that ofNa + . This is illustrated in Fig. 10 that shows model 
H + trajectories during a similar IMF turning. Here again, 
it can be seen that, for the test protons reaching the equator 
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Fig. 1L (Top) Model normalised density and (bottom) average energy of Na + ions in the noon-midnight meridian plane during turning of 
the IMF (from Bx = 0 to Bx — 20 nT in 20 s). The ions are launched from the high-latitude dayside sector (see text). The panels from left 
to right show snapshots at distinct limes of the magnetic transition. Black and grey lines show the magnetic field lines at the corresponding 
time and in the initial configuration, respectively. Density and energy are coded according to the color scales at right. 


before onset (trajectory coded in blue), similar results (e.g„ 
similar trajectory apex) are obtained regardless of the hemi- 
sphere of origin* In contrast, for protons located away from 
equator at onset, the trajectories significantly differ depend- 
ing upon the hemisphere of origin. Still, because of the 
much smaller gyroperiods of protons as compared to those of 
Na + , no significant ti change can be seen during IMF turn- 
ing (right panels of Fig. 10), which contrasts with the results 
in Fig. 8. 

To provide further insights into the response of ions to IMF 
turning, systematic trajectory computations were performed, 
considering test Na + launched from different altitudes (up to 
100km) and latitudes (in the 50° -90° range) in both hemi- 
spheres of the dayside sector. The test Na^ were initialized 
with 30° pitch angle and an energy of 10 eV comparable to 
that of exospheric neutrals (note that, due to the rapid cen- 
trifugal acceleration of the ions, little dependence was found 
upon these latter parameters). The test Na + were tracked 
according to a three-dimensional magnetospheric grid fea- 
turing bins of 0.05 Rm edge. Once the trajectory calcula- 
tions in the initial steady state (IMF Bx = 0) were performed, 
each bin was assigned a test particle with specific weight that 
was traced throughout turning of the IMF (from Bx = 0 to 
Bx — 20 nT in 20 s). The density (normalized to the maxi- 
mum value) and average energy in a given bin were subse- 
quently calculated by summation over the various lest ions 
occurring in that bin. The results of these computations are 
presented m Fig. 1 1 that shows the color-coded N density 
and energy at distinct times of the reconfiguration process. 

It can be seen in the leftmost panels of Fig. \ 1 that the 
patterns in the initial state are similar to those of Fig, 6, 
w r ith symmetrical flows above the polar cap and gradual cen- 


trifugal acceleration up to the keV range before entering the 
nightside plasma sheet. Most notably, moving from left to 
right in the bottom panels of Fig. 11, it is apparent that, in 
the course of the magnetic transition, the average energy of 
the Na^ population in the magnetospberic lobes rapidly in- 
creases up to several keVs. As discussed above (see Fig. 8), 
the present energization occurs in a nonadiabatic manner as a 
result of resonance between the induced electric field and the 
ion gyromotion. Because this short-lived nonadiabatic ener- 
gization depends upon gyroperiod versus held variation time 
scale, it preferentially affects heavy' ions (e.g., Na + ) and is 
negligible for protons. 

In the different panels of Fig, 1 1 * some ion flows with en- 
ergies of several keVs can also be seen north and south of 
the main Na^ population. Such flows result from peculiar 
ion drift paths as is illustrated in Fig. 12. It is apparent from 
this figure that, after transport into the nightside plasma sheet 
over the polar cap, the test Na + originating from the polar 
exosphere travels back into the dayside sector. It first mir- 
rors at high latitude near dusk, then in the cusp region near 
noon at which point it is convected back over the pole into 
the lobe with an energy of several keVs. A similar circula- 
tion scheme w p as invoked in the study of Seki et al. (1998) 
(see, e.g., their Fig, 6) to explain the Geotail observations 
of energetic (several keVs) G~ ions at large distances in the 
magnetospberic lobes. In addition, it is of interest to note 
that, at Mercury, the de-trapping zone of equatorial parti- 
cles (leading to the so-called “Shabansky branching”) ex- 
tends over most of the dayside sector (shaded area in the 
bottom left panel of Fig. 12)* That is, because of the exis- 
tence of a field minimum at high latitude in the cusp region, 
the mirror force points away from equator near the frontside 
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Fig, 12* Model trajectory from polar exosphere to high- 
latitude magnetospheric lobe via nights ide plasma sheet: (left) tra- 
jectory projections in tire noon-midnight meridian plane and in the 
equatorial plane, (right) energy and magnetic moment (normalized 
to the initial value) as a function of time. In the bottom left panel, 
the thick Line and shaded area show the magnetopause and the re* 
gion of cuspward mirror force, respectively. 


magnetopause, instead of equatorward elsewhere. Accord- 
ingly, (nearly) equatorially trapped particles are released to- 
ward high latitudes when they approach the frontside mag- 
netopause, which provides a mechanism for populating the 
high- latitude region with energetic particles (Delcourt and 
Sauvaud, 1999). While at Earth, this de- trapping zone is lim- 
ited to the very vicinity of the frontside magnetopause {see, 
e.g.. Fig. 1 of Delcourt and Sauvaud, 1999), it extends over 
a wide domain at Mercury due to the small planetary field, 
hence hampering stable trapping at equator through most of 
the dayside magnetosphere. 


5 Conclusions 

The three-dimensional single-particle computations per- 
formed demonstrate that the IMF orientation has a promi- 
nent impact on the large-scale ion circulation at Mercury. A 
nonzero IMF Bx should lead to asymmetry between northern 
and southern magnetospheric lobes with a more pronounced 
E x B related centrifugal acceleration in one of them; hence, 
a more rapid filling of this very lobe and distinct accessible 
regions in the magnetotail depending upon the hemisphere 
of origin of the ions. The simulations also reveal that, in a 
similar manner to compression events (e.g + , Delcourt et a T, 
2010), IMF turning at Mercury should lead to resonant heat- 
ing of the ions under the effect of the induced electric field. 
That is, for ions that have gyration periods comparable to the 
field variation time scale, a possibly significant energization 


(up to the keV range and above) is obtained, characterized 
by prominent magnetic moment increase. Such a short-lived 
nonadiabatic heating of ions at high latitudes appears spe- 
cific to the hermean environment that is characterized by a 
weak magnetic field as well as rapid field variations. During 
IMF turning, it should lead to hemisphere-dependent loading 
of the magnetospheric lobes with heavy energetic material, 
a feature that may be elucidated with the forthcoming MES- 
SENGER data. 
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